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FIELD 1 SYMBOL TABLE STORAGE FOR PALD (DEC-D8-ASAA-LA) 


DECUS Program Library Write-up DECUS No. 8-376 


PALD provides for storing symbols on the system device, when 
the number of symbols exceeds 1o0 _ When using DECtape, assembly and 
listing time can be excessive, due to the delays involved in the symbol 
table search on the mag. tape. 

If Field 1 is available, the DECtape delays can be eliminated 
by storing the symbol table in Field 1, reducing assembly and listing 
time considerably. The saving in time will depend on the number of 
external symbols which are generated by the program. For a program 
requiring three blocks of external symbols, the saving in assembly time 
(compared to a DECtape symbol table) was one half. The larger the symbol 
table, the more time will be saved. 'With a disc as system device, the 
time saving is not significant. 

Twenty blocks (of 200 Q words) are reserved in Field 1 for the 
symbol table, replacing . SYM on the system device. 

Later versions of PALD may require some changes to the overlay. 
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STATPAC REVISIONS FOR PDP-S/l AND TSS/8 


DECUS Program Library Write-up 


DECUS NO. 8-406 


ABSTRACT 

This package contains 11 programs from the original PDP-10 Dartmouth BASIC Statistical 
Package which have been revised for the PDP-8/1 and TSS/8. The documentation consists 
of a description and listing of each of the programs. 








10 ’ BERKSHIRE COMMUNITY COLLEGE USER LIBRARY 

20 ’ MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 
30 ’ DARTMOUTH PROGRAM NAME: 3TAT01*** 

40 ’ MODIFIED: JULY 2 1970 
90 * 

100* B.C.C. LIBRARY PROGRAM NAME: STAT01 

110* 

120' DESCRIPTION--COMPUTES THE MEAN, VARIANCE, STANDARD DEVIATION, 
130* AND THE STANDARD ERROR OF THE MEAN FOR ONE OR MORE 
140’ SETS OF DATA. 

150* 

160' SOURCE--UNKNOWN 

170* 

180' INSTRUCTIONS - - THE DATA START IN LINE NUMBER 900 AND CONTINUE 
190* NOT PAST LINE 909. THE NUMBER 999999 MUST BE PLACED AT 
200‘ THE END OF THE LAST DATUM OF EACH SET. SAMPLE DATUM ARE IN LINE 
210’ 900. 

220 * 

230' 

240* ******* MAIN PROGRAM ********** 

250’ 

260 READ X 

270 IF X = 999999 THEN 320 

280 LET S = S + X 

290 LET SI = SI + X*X 

300 LET N = N + 1 

310 GO TO 260 

320 LET M = S/N 

330 LET V = (N*S1 -S*S)/N/CN-1) 

340 LET S3 = SQRCV) 

350 LET S4 = SQRCV/N) 

360 PRINT"NUMBER", "SUM SUM-0F-SQUARES","MEAN","VARIANCE" 

370 PRINT N, S, SI, M, V 

380 PRINT 

390 LET Q = S3/M 

400 PRINT "STD. DEV.", "STD. ERROR OF THE MEAN","COEFF. OF VARIATION" 

410 PRINT S3, S4, " ", Q 

420 PRINT 

430 PRINT 

440 READ X 

450 IF X = 999999 THEN 500 

460 LET S=0 

470 LET S1=0 

480 LET N=0 

490 GOTO 280 

500 STOP 

900 DATA 1,2,3,4,5,6,7,8,9, 10,999999,4,3,9, 15,21,6,999999 
910 DATA 999999 
920 END 

READY 
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10 ' BERKSHIRE COMMUNITY COLLEGE USER LIBRARY 

20 ’ MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 
30 ’ DARTMOUTH PROGRAM NAME: STAT02*** 

40 ' MODIFIED: JULY 2 1970 
90 • 

100 ' B.C.C LIBRARY PROGRAM NAME: STAT02 

110 * 

120’ DESCRIPTION--COMPUTES THE MEANS* VARIANCES* AND T-RATIO 

130* FOR TWO GROUPS OF DATA*UNPAIRED. THIS PROGRAM ASSUMES THAT THE 

140* TWO GROUPS HAVE THE SAME VARIANCE. 

150* 

160* SOURCE--THOMAS E. KURTZ* KIEWIT COMPUTATION CENTER 

170* 

180* INSTRUCTIONS--PUT DATA IN LINES 900 THROUGH 998. 

190’ END THE FIRST SERIES OF DATA WITH 999999 AND THEN TYPE 
200' THE SECOND SERIES* AGAIN ENDING WITH 999999. 

210* SAMPLE DATA ARE IN LINES 900 AND 910. 

220 * 

230* 

240* 

250* ****** MAIN PROGRAM ********** 

260 DIM SC2)*ZC2)*NC2)*MC2)*VC2)*DC2) 

270 LET P = 1 
280 LET S = 0 
290 LET S2 = 0 
300 LET N = 0 
310 READ X 

320 IF X = 999999 THEN 370 

330 LET S = S + X 

340 LET S2 = S2 + X*X 

350 LET N = N + 1 

360 GO TO 310 

370 LET S(P) = S 

380 LET Z(P> = S2 

390 LET NCP) = N 

400 IF P = 2 THEN 430 

410 LET P = 2 

420 GO TO 280 

430 REM NOW WE PRINT THE ANSWERS 

440 PRINT "GROUP”* "NUMBER"* "MEAN"* "VARIANCE"* "STD. DEV." 

450 PRINT 

460 FOR I = 1 TO 2 

470 LET MCI) = SC D/NCI) 

480 LET VCI) = CNCI)*ZCI) - SC I)*SCI))/NCI)/CNCI) - 1) 

490 LET DC I) = SQRCVCI)) 

500 PRINT I* NCI)* MCI)* VCI)* DC I) 

510 NEXT I 

520 LET 0 = VC 1)*CNC1)-l) + VC2)*CNC2)-l) 

530 LET Q = Q/C NCI) + N C 2) -2 ) * C 1/NCI) + 1/NC2) ) 

540 LET R = MCI) -MC2) 

550 LET W = SORCQ) 

560 PRINT 

570 PRINT "MEAN DIFF. "* "VAR. DIFF. "* "STD. DEV. DIFF." 

58U PRINT R* Q* W 
590 PRINT 

600 PRINT "T RATIO"* R/W* "ON"; NCl)+NC2)-2* "DEGREES OF FREEDOM." 
610 STOP 

900 DATA 160* 160* 140* 190*999999 

910 DATA 1 17* 145* 147* 120* 150* 120*999999 

999 END 


READY 


3 







B^V^§iH s 5^Leiig^AVI§^l;^ program 

DARTMOUTH PROGRAM NAME: STAT03*** 

MODIFIED: JULY 2 1970 


DESCRIPTION--COMPUTES THE MEANS# VARIANCES# AND T-RATI3 
FOR TWO GROUPS OF UNPAIRED DATA. THIS PROGRAM ASSUMES THAT 
GROUPS MAY HAVE UNEQUAL VARIANCES. 


INSTRUCTIONS--PUT DATA IN LINE 900 AND FOLLOWING. 
MAKE SURE THE DATA LINE NUMBERS DO NOT EXCEED 998. 
END THE FIRST SERIES OF DATA WITH 999999# AND THEN 
TYPE IN THE SECOND SERIES# AGAIN ENDING WITH 999999. 
SAMPLE DATA ARE IN LINES 900 AND 910. 


* * 


* 


THEN 370 


X*X 


******* MAIN PROGRAM 
DIM SC2)#ZC2)#NC2)#MC2)# VC2)#DC2) 

LET P = 1 
LET S = 0 
LET S2 = 0 
LET N = 0 
READ X 

IF X = 999999 
LET S = S + X 
LET S2 = S2 + 

LET N = N + 1 
GO TO 310 
LET SCP) = S 
LET ZCP) = S2 
LET NCP) = N 
IF P = 2 THEN 
LET P = 2 
GO TO 280 

REM NOW WE PRINT THE ANSWERS 
PRINT "GROUP"# "NUMBER"# "MEAN"# 

PRINT 

FOR I = 1 TO 2 
LET MCI) = SC D/NCI) 

= CN CI)*Z Cl) - SCI)*SCI))/NCI)/CNCI) - 
= SQRCVCI)) 

NCI)# MCI)# VCI)# 


* 


£8 

30 

40 

90 ’ 

100 * B.C.C. LIBRARY PROGRAM NAME: STAT03 
110 * 

120 ' 

130’ 

140* 

150' 

160' SOURCE--THOMAS E-. KURTZ# KIEWIT COMPUTATION CENTER 
170* 

180' 

190' 

200 ' 

210 ’ 

220 ' 

230' 

240* 

250' 

260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
900 
910 
999 


THE 


* 


430 


•VARIANCE' 


LET VCI) 
LET DC I) 
PRINT I# 
NEXT I 


STD. DEV. 


1 ) 


Q 

W 

R 

K 

D 


LET 
LET 
LET 
LET 
LET 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
STOP 
DATA 
DATA 
END 


VCI) / 
SQRCQ) 
MCI) - 
VCI) / 


DC I) 


NCI) + VC 2) / NC 2) 


MC 2) 
NCI) 


= 1/C K*K/CNC1) 


Q 

- 1 ) 


+ C1-K)*C1-K)/CNC2) - 1) ) 


"MEAN 

R# Q# 


DIFF. 

W 


•VAR. DIFF."# "STD. DEV. DIFF." 


"T RATIO"# R/W # 


OF FREEDOM. 


1 60# 
117 # 


1 60# 
145# 


1 40# 
1 47# 


190# 

120 # 


999999 
150# 1 20j 


999999 


** 


READY 
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10 * BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 ' MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 • DARTMOUTH PROGRAM NAME: S TAT04*** 

40 • MODIFIED: JULY 2 1970 
90 * 

100 ' B.C.C. LIBRARY PROGRAM NAME: STAT04 

110* 

120* DESCRIPTION--COMPUTES CHI SQUARE STATISTICS FOR 
130* 2 BY 2 CONTINGENCY TABLES. 

140* 

150* SOURCE--THOMAS E. KURTZ* KIEWIT COMPUTATION CENTER 

1 60* 

170’ INSTRUCTIONS--DATA IS PUT IN LINE 900 AND FOLLOWING. 

180* ENTER THE TABLES IN NORMAL ORDER WITH FOUR ENTRIES OF 
190* TABLE 1 FIRST* THEN THE 4 ENTRIES FOR TABLE 2* ETC. 

200* DATA LINE NUMBERS SHOULD NOT EXCEED 998. SAMPLE DATA ARE 
210* IN LINE 900. 

220 * 

230* 

240' * * He * * * MAIN PROGRAM ********** 

250' 

260 PRINT "TABLE"* " "* "CHI SQUARE" 

270 PRINT 

280 READ A* B* C* D 

290 LET N=A+B+C+D 

300 LET E = A*D - B*C 

310 LET E=ABS< E) 

320 LET F = <E - N / 2 ) 

330 LET G = N * F * F 

340 LET R1 = A + B 

350 LET R2 = N - R1 

360 LET Cl = A + C 

370 LET C2 = N - Cl 

380 LET X = G / R1 / R2 / Cl / C2 

390 PRINT A* B 

400 PRINT C* D* X 

410 GO TO 270 

900 DATA 6* 0* 2* 6* 6* 0* 1* 3 

999 END 

READY 


5 







10 * BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 * MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 ' DARTMOUTH PROGRAM NAME: STAT08*** 

40 * MODIFIED: JULY 2 1970 
90 ' 

100 ' B» C» C • LIBRARY PROGRAM NAME: STAT08 
110 ' 

120* DESCRIPTION--C0MPARES TWO GROUPS OF DATA USING THE MEDIAN TEST. 

1 30* 

140’ SOURCE--THOMAS E. KURTZ, KIEWIT COMPUTATION CENTER 
150' 

1 60' INSTRUCTIONS--PUT DATA IN LINE 900 AND FOLLOWING. 

170' THE FIRST DATUM IS THE NUMBER OF ENTRIES IN THE FIRST 

180* GROUP, THEN ENTER THE NUMBER OF ENTRIES IN THE SECOND GROUP, 

190* THEN THE FIRST GROUP ITSELF IS ENTERED,AND THEN THE SECOND GROUP. 
200* THE PROGRAM PRINTS OUT THE CHI SQUARE STATISTIC OF A 2 BY 2 TABLE 
210' ON 1 DEGREE OF FREEDOM. DATA SHOULD NOT EXCEED LINE 998. 

220’ SAMPLE DATA ARE IN LINE 900. 


230* 

240' 

250' ****** MAIN PROGRAM 

260' 

265 DIM AC 100) 

270 READ M, N 

280 LET Ml = M + N 

290 LET M2 = INTC Ml/2 ) 

295 FOR Zl= 1 TO Ml 
300 READ ACZ1) 

305 NEXT Z1 

310 LET L = 1 

320 LET U = M 

330 GOSUB 630 

340 LET L = M + 1 

350 LET U = Ml 

360 GOSUB 630 

365 LET X = 0 

370 LET Y = 0 

380 LET I = 1 

390 LET J = M + 1 

400 FOR K = 1 TO M2 

410 IF AC I ) < AC J) THEN 470 

420 LET Y = Y + 1 

430 LET J = J + 1 

440 IF J <= Ml THEN 520 

450 LET X = M2 - Y 

460 GO TO 530 

470 LET X = X + 1 

480 LET 1=1+1 

490 IF I <= M THEN 520 

500 LET Y = M2 - X 

510 GO TO 530 


520 
530 
540 
550 
560 
57 0 
580 
590 
600 
610 
620 
630 


NEXT K 
LET U = 
LET V = 
LET Z = 
LET T = 
LET C2 
PRINT 
PRINT 
PRINT 
PRINT 
STOP 
FOR I = 


M - X 
N - Y 
X*V - Y*U 

M1 *C ABSCZ) - Ml/2 ) t 2 
T/M/N/M2/CM1- M2) 
"TWO SAMPLE MEDIAN TEST." 
"GROUP 1 "5 x; U 

"group 2 y; v 

"CHI-SQUARE = "i C2 


6 


* * 


* 


* * * * 


* * * 


1 TO U 


L 





640 LET X = ACL) 

650 LET Q = L 

660 FOR J = L + 1 T0 U • I + 1 

670 IF X >= ACJ) THEN 700 

680 LET X = ACJ) 

690 LET Q = J 

700 NEXT J 

705 LET J = U-I+l ’THIS STATEMENT PUT IN TSS/8 BECAUSE 0F BUG IN BASIC 
710 IF J = Q THEN 740 

720 LET ACQ) = ACJ) 

730 LET ACJ) = X 

740 NEXT I 
750 RETURN 

900 DATA 4, 6, 160, 160, 140, 190, 117, 145, 147, 120, 150, 120 

999 END 

READY 
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10 • BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 * MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 ' DARTMOUTH PROGRAM NAME: STAT9X*** 

40 • MODIFIED: JULY 2 1970 
90 * 

100 ’ B.C.C. LIBRARY PROGRAM NAME: STAT9X 

110* 

120' DESCRIPTION-“COMPUTES THE SLOPE AND OTHER STATISTICS FOR A 
130* LINEAR REGRESSION WITH SEVERAL Y VALUES FOR EACH X VALUE. 

1 40' 

150* S0URCE--THOMAS E» KURTZ* KIEWIT COMPUTATION CENTER 

1 60* 

170' INSTRUCTIONS--PUT DATA STARTING IN LINE 900. 

180' FIRST DATA IS K, THE NUMBER OF DIFFERENT X-VALUES. THEN 
190' FOR EACH OF THE K GROUPS-ENTER THE NUMBER IN THAT GROUP* 
200* THEN THE COMMON X-VALUE* THEN THE Y-VALUES FOR THAT GROUP. 
210* SAMPLE DATA ARE IN LINES 900-950. 

220 * 

230* 

240* ****** MAIN PROCRAM * * * * * * * * * 

250 READ K 

260 FOR J = 1 TO K 

270 READ NO 

280 LET N=N+NO 

290 READ X 

300 FOR I = 1 TO NO 

310 READ Y 

320 LET XI = XI + X 

330 LET Y1 = Y1 + Y 

340 LET X2 = X2 + X*X 

350 LET Y2 = Y2 + Y*Y 

360 LET Z = Z + X*Y 

370 NEXT I 

380 NEXT J 

390 LET SI = N*X2 - X1*X1 

400 LET S2 = N*Z - Xl*Yl 

410 LET B = S2/S1 

420 LET Y3 = Yl/N 

430 LET X3 = Xl/N 

440 LET B1 = Y3 - B*X3 

450 LET N1 = N - 1 


460 LET N2 = N1 - 1 

470 LET S3 = CY2 - Y1*Y3 - B*S2/N> 

480 LET S4 = S3/N2 

490 PRINT "NUMBER = "J N* "SLOPE = "*B 
500 PRINT "MEAN OF X = "J Xl/N* " OF Y 
510 PRINT "Y-INTERCEPT = "B1 
520 PRINT "SUM-OF-SQUARES"* "TOTAL"*Y2 
530 PRINT " "* "MEAN "* Y3*Y1 


"i Yl/N 


540 PRINT " 
550 PRINT " 
560 PRINT " 
570 PRINT " 
580 PRINT " 
590 PRINT " 
600 PRINT " 
610 PRINT " 
620 PRINT " 
630 PRINT " 
640 STOP 
900 DATA 5 
910 DATA 3* 


"* "SLOPE "* B*S2/N 
"* "RESIDUAL "* S3 
STANDARD DEVIATIONS" 

", "X "* SQRCS1/N/N1) 

", "Y "* SQRC <Y2-Y1*Y3)/N1) 
"* "ERROR "* SQRCS4) 

"* "Y-BAR "* SQRC54/N) 

"* "SLOPE "* SQRC S4/S1*N) 

"* "Y-INTER."* SQRCS4*( 1 /N + 
F-RATIO FOR SLOPE = " B*S2/N/34 


4* 6* 8* 7 


920 DATA 2* 7* 1* 3 


X3*X3/S1*N) 


) 




930 

940 

DATA 

DATA 

5, 

It 

6t 

It 

It 

2 

2t 

3t 

At 

5 

950 

999 

DATA 

END 

At 

5t 

21 

At 

6 1 

3t 

10 


READY 







10 • BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 ’ MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 * DARTMOUTH PROGRAM NAME: STATU*** 

40 * MODIFIED: JULY 2 1970 
90 * 

100 ’ B.C.C. LIBRARY PROGRAM NAME: STAT11 
110’ 

120' DESCRIPTION--COMPUTES THE SPEARMAN RANK CORRELATION COEFFICIENT 
130’ FOR TWO SERIES OF DATA. 

140’ 

150* SOURCE--THOMAS E. KURTZ* KIEWIT COMPUTATION CENTER 
1 60’ 

170* INSTRUCTIONS--PUT DATA IN LINE 900 AND FOLLOWING. 

180* ENTER DATA IN THE FOLLOWING ORDER: FIRST THE COMMON LENGTH 
190’ OF EACH SERIES* THEN THE ELEMENTS OF THE SERIES BY GROUPS 
200’ THAT IS X1*Y1*X2*Y2* .. ..ETC. 

210* SAMPLE DATA ARE IN LINE 900. 

220 ’ 

230’ 

240’ ******* MAIN PROGRAM ******** 

250* 

260 DIM AC 100)* BC100) 

270 READ N 

280 FOR I = 1 TO N 

290 READ AC I)* BCI) 

300 NEXT I 

310 FOR I = 1 TO N - 1 
320 LET X = AC 1) 

330 LET Y = BC1) 

340 LET Q = 1 

350 FOR J = 2 TO N - I + 1 

360 IF X >= ACJ) THEN 400 

370 LET X = AC J) 

380 LET Y = BCJ) 

390 LET Q = J 

400 NEXT J 

410 LET J = N-I+l 

411 IF 9=J THEN 460 

420 LET-ACQ) = ACJ) 

430 LET ACJ) = X 

440 LET BCQ) = BCJ) 

450 LET BCJ) = Y 

460 NEXT I 

470 REM NOW WE REPLACE THE ELEMENT BY ITS AVERAGE RANK 
480 FOR I = 1 TO N 
490 LET X = AC I) 

500 FOR J = I + 1 TO N 

510 IF X <> ACJ) THEN 540 

520 NEXT J 

530 LET J = N + 1 

540 LET Z = CJ + I - l)/2 

550 FOR K = I TO J - 1 

560 LET ACK) = Z 

570 NEXT K 

580 LET I = J - 1 

590 NEXT I 
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600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

900 

910 


IF P = 1 THEN 690 

REM NOW WE INTERCHANGE A AND B AND SORT AGAIN. 

FOR I = 1 TO N 
LET T = AC I) 

LET AC I) = BCI) 

LET BCI) = T 
NEXT I 
LET P = 1 
GO TO 310 

REM NOW WE START COMPUTING THE CORRELATION COEFFICIENT. 
FOR I = 1 TO N 
LET SI = SI 


LET S2 = S2 - 
LET S3 = S3 J 
NEXT I 

LET S = N * N * 

LET D1 = N * SI 
LET D2 = N * S2 
LET D = D1 * D2 
IF D <> 0 THEN 820 


AC I ) 
BCI) 
AC I ) 

CN+1 ) 

- S 

- S 


AC I ) 
BCI) 
BCI) 


* CN+1) / 4 


LET R 
GO TO 
LET R 
PRINT 
PRINT 
STOP 
DATA 
END 


= 0 
830 

= CN * S3 - S)/SQRC D) 
"SPEARMAN RANK CORRELATION 


COEFFICIENT" 


5* 430, 56, 500, 61, 520, 78, 540, 71, 56, 82 


READY 
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10 ' BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 * MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 * DARTMOUTH PROGRAM NAMES STAT20*** 

40 * MODIFIED: JULY 2 1970 
90 ' 

100 * B.C.C. LIBRARY PROGRAM NAME: STAT20 
110 REM JUNE 11, 1970 

120 REM DESCRIPTION-- MULTIPLE LINEAR REGRESSION ACCORDING TO 
130 REM EFROYSON'S ALGORITHM 
140 REM 

150 REM SOURCE-- RALSTON AND WILF,"MATHEMATICAL METHODS FOR DIGITAL 
160 REM COMPUTERS’ 1 , P. 191 

165 REM PROGRAMMED BY THOMAS E. KURTZ, KIEWIT COMPUTATION CENTER 
170 REM 

180 REM INSTRUCTIONS--ENTER DATA STARTING IN LINE 1610. 

190 REM FIRST DATA IS N, THE NUMBER OF INDEPENDENT VARIABLES, THEN 
200 REM P, THE NUMBER OF DEPENDENT VARIABLES, THEN M, THE NUMBER 
210 REM OF DATA SETS, THEN FI, THE VALUE OF F FOR ENTERING 
220 REM A VARIABLE, AND THEN F2, THE VALUE OF F FOR REMOVING A 
225 REM VARIABLE. THEN ENTER THE DATA BY SETS, INDEPENDENT VARIABLES 
230 REM BEFORE DEPENDENT VARIABLES. IF M>50 OR P>7 OR N>8 THE 
240 REM DIM STATMENTS IN LINE 300 SHOULD BE CHANGED. SAMPLE DATA 
250 REM ARE IN LINES 1610-1760. 

260 REM BE SURE TO REMOVE SAMPLE DATA BEFORE RUNNING PROGRAM. 

270 REM 

280 REM ******* MAIN PROGRAM ******* 

290 REM 

300 DIM S( 6, 6) 

301 DIM A(20, 4) 

302 DIM MC8),DC8),BC7),EC7) 

310 READ N, P, M, FI, F2 
320 LET T1 = IE-6 

330 LET N1 = N + P 
340 LET N2 = N1 + N 
350 LET K = 1 
360 LET D1 = M - 1 

370 REM MAT READ ACM,N1> 

371 FOR 19 = 1 TO M 

372 FOR 18 = 1 TO N1 

373 READ ACI9,18) 

375 NEXT 18 

376 NEXT 19 

380 FOR I = 1 TO M 
390 LET ACI,0)=1 

400 NEXT I 
410 PRINT "MEANS " 

420 FOR I = 0 TO N1 
430 FOR J = I TO N1 

440 LET S = 0 

450 FOR L = 1 TO M 

460 LET S = S + ACL,I) * ACL,J> 


470 NEXT L 

480 LET SCI,J) = S 

490 NEXT J 

500 LET MCI) = SCO,I) / SC0,0> 

510 IF I = 0 THEN 530 

520 PRINT MCI), 

530 NEXT I 
540 PRINT 

550 PRINT 

560 PRINT "STANDARD DEVIATIONS" 

570 LET Ml = M * D1 
580 FOR I = 1 TO N1 
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C M*SCI#J) 


SCO# I )*SCO# J) )/Ml 





590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

7 60 
770 
780 
790 
800 
810 
820 
830 
840 
850 

8 60 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 

1000 

1010 

1020 

1030 

1031 

1032 

1033 
1040 
1050 
1060 
1070 
1030 
1090 
1100 
1110 
1120 
1130 
11 40 
1150 
1 1 60 
1170 
1180 


FOR J = I TO N 1 
LET SCI#J> = 

NEXT J 

LET DC I) = SQRC SC I#I) ) 

PRINT DC I) # 

NEXT I 
PRINT 
PRINT 

PRINT "CORRELATION COEFFICIENTS" 

FOR I = 1 TO N2 
FOR J = I TO N2 

IF J > N 1 THEN 740 

LET SCI#J) = SCI#J) / DC I) / DC J) 
LET SCJ#I> = SCI#J> 

GO TO 800 

IF I <> J - N1 THEN 780 
LET SCI#J) = 1 
LET SCJ#I> = -1 
GO TO 800 

LET SCI#J> = 0 
LET SCJ#D = 0 
NEXT J 
NEXT I 

FOR I = 1 TO N1 

FOR J = 1 TO N1 

PRINT SCI#J)# 

NEXT J 
PRINT 
PRINT 
NEXT I 
PRINT 

LET K1 = K + N 

PRINT "DEPENDENT VARIABLE"#K 
FOR I = 1 TO N 
LET BCI) = 0 
NEXT I 

LET S8 = DCK1) * SQRC SCK1#K1> / D1 ) 
LET VI = 1E35 
LET V2 = 0 
LET N3 = 0 
LET N4 = 0 
FOR I = 1 TO N 

IF ABSC SC I#I) ) <= T1 THEN 1140 
LET VO = SC I# K1) * SCK1#I) / SCI#I) 
LET E9= SGN C VO) + 2 
IF E9 = 1 THEN 1080 
IF E9 = 2 THEN 1 140 
IF E9 = 3 THEN 1040 

IF VO <= V2 THEN 1140 
LET V2 = VO 
LET N4 = I 
GO TO 1140 
LET II = I + N1 

LET BCD = SC I 1 # K 1 ) * DCK 1) / DC I) 
LET EC I) = S8 / DC I) * SQRC SCI1#I1) 
IF ABSC VO ) >= ABSC VI ) THEN 1140 
LET VI = VO 
LET N3 = I 
NEXT I 
LET S = 0 
FOR I = 1 TO N 

LET S = S + BCD * MCI) 

NEXT I 


) 
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1190 LET BCO) = M C K 1) - S 

1200 PRINT "INDEX ", "B ", "STD. DEV.", "T-RATIG 
1210 PRINT 0, BCO) 

1220 F0R I = 1 T0 N 

1230 IF BCI) = 0 THEN 1250 

1240 PRINT I, BCI), EC I), BCI) / EC I) 

1250 NEXT I 

1260 PRINT "STANDARD ERROR 0F Y ", S8 * SQRCM-1) 

1270 PRINT 

1280 LET F = ABSCV1) * D1 / SCK1,K1) 

1290 IF F < F2 THEN 1360 

1300 LET F » V2 * (D1 - 1) / C SCK1,K1) - V2 ) 

1310 IF F <= FI THEN 1570 
1320 LET Q = N4 
1330 LET D1 = D1 — 1 

1340 PRINT "VARIABLE ENTERING ", Q 

1350 G0 T0 1390 
1360 LET Q = N3 
1370 LET D1 = D1 + 1 

1380 PRINT "VARIABLE LEAVING ", Q 

1390 PRINT "F-LEVEL ", F 

1400 LET Y = 1 / SCQ,Q) 

1410 F0R J = 1 T0 N2 

1420 LET SCQ,J) = SCQ,J) * Y 

1430 NEXT J 

1440 LET Y = -Y 

1450 F0R I = 1 T0 N2 

1460 IF I = Q THEN 1540 

1470 LET X = -SC I,Q) 

1480 F0R J = 1 T0 N2 

1490 IF J = Q THEN 1520 

1500 LET SCI,J) = SCI,J) + X * SCQ,J) 

1510 G0 T0 1530 

1520 LET SCI,J) = SCI,J) * Y 

1530 NEXT J 

1540 NEXT I 

1550 LET SCO,Q) = -Y 

1560 G0T0 950 

1570 LET K = K + 1 

1580 IF K <= P THEN 900 

1590 PRINT " ***** END OF PROBLEM *****•• 

1600 ST0P 

1610 DATA 3, 1,9, 2. 5, 2. 5 
1620 DATA 24769,34181,35516,29566 
1630 DATA 34181,35516,29566,37009 
1640 DATA 35516,29566,37009,34542 
1650 DATA 29 
566,37009,34542,41485 


1 660 

DATA 

37009,34542,41485, 

49972 

1670 

DATA 

34542,41485,49972, 

53944 

1680 

DATA 

41485,49972,53944, 

47062 

1 690 

DATA 

49972,53944,47062, 

47284 

1700 

DATA 

53944,470 62,47284, 

57765 

1710 

DATA 

1, 1, 1, 1 


1720 

DATA 

3,4 


1730 

DATA 

2,2,1,0 


1740 

DATA 

1,3,4 


1750 

DATA 

3,3,0,1 


1755 

DATA 

1,2,3,4 


17 60 

DATA 

o 

* 

o 

* 

* 

sr 


1999 

END 




READY 


14 






STAT21 


MULTIPLE LINEAR REGRESSION - FANCY OUTPUT 


DESCRIPTION - Computes one or more multiple linear regressions on a 
batch of data. 

INSTRUCTION - Place data in line 1940 and following. First data is 
N (number of data sets or observations), then V (num¬ 
ber of variables in data batch), then G (number of re¬ 
gressions to be performed). Then enter the data by the 
data set, that is enter first the values for all vari¬ 
ables at observation 1, then the values for variables 
at obvervation 2, etc. Variables may be entered in any 
order, but the order must be the same for each data set. 
The position of a variable in entering the data is the 
index of the variables. Thus, if there are four vari¬ 
ables entered in the order X, Y, Z, W, the index of Z 
is 3. Next, for each regression enter (in order) H 
(the number of the regression), then K (total number of 
independent variables in the regression), then PI 
(enter 1 if you want the residuals printed out, other¬ 
wise enter 0), P2 (enter 1 of you want the residuals 
printed out, otherwise enter 0)., and then enter the in¬ 
dexes of the independent variables followed by the 
index of the dependent variable. If N > 9 or V > 4, 
then the DIM statements in line 430 and 440 must be 
changed. Sample data are in lines 1940 through 1992. 

Be sure to remove the sample data before running the 
program. 
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10 ' BERKSHIKE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 • MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 ' DARTMOUTH PROGRAM NAME: STAT21*** 

40 * MODIFIED: JULY 2 1970 
90 * 

100 • B.C.C. LIBRARY PROGRAM NAME: STAT21 
110 REM 5/27/70 
120 REM 

130 REM DESCRIPTION-- COMPUTES ONE OR MORE MULTIPLE LINEAR 
140 REM REGRESSIONS ON A BATCH OF DATA. 

145 REM 

150 REM SOURCE--REVISED 11/16/67 BY PROF. GERALD CHILDS. 

160 REM 

170 REM INSTRUCTIONS--PLACE DATA IN LINE 1960 AND FOLLOWING. 

180 REM FIRST DATA IS N (NUMBER OF DATA SETS OR OBSERVATIONS)* 

190 REM THEN V (NUMBER OF VARIABLES IN DATA BATCH)* THEN G 

200 REM (NUMBER OF REGRESSIONS TO BE PERFORMED)* THEN ENTER THE 

210 REM DATA BY THE DATA SET* THAT IS ENTER FIRST THE VALUES FOR 

220 REM ALL VARIABLES AT OBSERVATION 1, THEN THE VALUES FOR THE 

230 REM VARIABLES FOR OBSERVATION 2* ETC. VARIABLES MAY BE 

240 REM ENTERED IN ANY ORDER* BUT THE ORDER MUST BE THE SAME FOR EACH 

250 REM DATA SET. THE POSITION OF A VARIABLE IN ENTERING THE DATA IS 

260 REM THE INDEX OF THE VARIABLE. THUS* IF THERE ARE 4 VARIABLES 

270 REM ENTERED IN ORDER X*Y*Z*W THE INDEX OF Z IS 3. NEXT* 

280 REM FOR EACH REGRESSION ENTER (IN ORDER) H (THE NUMBER OF THE 
290 REM REGRESSION)* THEN K (TOTAL NUMBER OF INDEPENDENT VARIABLES 
420 REM ********* MAIN PROGRAM ************ 

430 DIM X(4*4)*A(4*4)*D(9*4)*Y(4)*M(4)*S(4) 

440 DIM T(4)*B(4)*U(4*4)*R(4*4)*C(4*4)*0(9)*E(4) 

450 READ N* V* G 

451 REM MAT READ D(N*V) 

452 FOR I = 1 TO N 
454 FOR J = 1 TO V 
456 READ D(I*J> 

458 NEXT J 

459 NEXT I 

470 FOR 1=1 TO N 
480 LET D( I * 0 ) = 1 
490 NEXT I 

500 READ H* K* PI* P2 
510 LET M = K + 1 

514 REM MAT READ E(M) 

515 FOR Z = 1 TO M 
520 READ E(Z) 

525 NEXT Z 

530 PRINT "**REGRESSI3N NUMBER"H": DEPENDENT VARIABLE IS"E(M) 

540 PRINT 

550 IF H>1 THEN 700 

560 FOR 1 = 0 TO V 

570 FOR J=0 TO V 

580 LET X=0 

590 FOR L= 1 TO N 

600 LET X=X+D(L*I)*D(L*J) 

610 NEXT L 
620 LET X(I*J)=X 
630 LET C(I*J)=X 

640 NEXT J 

650 LET T(I)=X(0*I)/X(0*0) 

660 LET B(I) = 0 
670 IF 1=0 THEN 690 

680 LET B( I)=S9R(X(I * I)/(N-1)-X(0*I)*X(0*I)/(N*(N-1))) 

690 NEXT I 

700 PRINT ’'INDEX"* "MEANS"* "STANDARD DEVIATIONS" 16 





710 FOR 1 = 1 TO M 
720 LET M(I) = T C E CI )) 

730 LET SCI)=BCE(I)) 

740 PRINT EC I ) , MCI), SCI) 

750 NEXT I 
760 PRINT 
770 PRINT 

780 PRINT "CORRELATION COEFFICIENTS" 

790 IF H>1 THEN 850 
800 FOR 1=1 TO V 
810 FOR J=1 TO V 

820 LET RCI,J)=CN*XCI,J)-XCO,I)*XC0,J))/CN*CN-1)*3CI)*BCJ) 

830 NEXT J 

840 NEXT I 

850 FOR 1=1 TO M 

8 60 FOR J=1 TO M 

870 LET UCI,J)=RCECI),ECJ)) 

880 PRINT UCI,J), 

890 NEXT J 

900 PRINT 

910 PRINT 

920 NEXT I 

930 PRINT 

940 LET ECO)=0 

950 FOR 1 = 0 TO K 

960 LET YCI)=CCECI),ECM)) 

970 FOR J=0 TO K 

980 LET XCI,J)=CCECI),ECJ)) 

990 NEXT J 

1000 NEXT I 
1010 FOR 1 = 0 TO K 
1020 FOR J=0 TO K 
1030 IF I<>J THEN 1060 
1040 LET AC I,J) = 1 
1050 GO TO 1070 
1060 LET AC I,J)=0 

1070 NEXT J 

1080 NEXT I 
1090 FOR 1 = 0 TO K 
1100 IF XCI,I)<lE-6 THEN 1840 
1110 LET YCI) = YCI)/XCI,I ) 

1120 FOR J=0 TO K 

1130 LET AC I,J)=ACI,J)/XCI,I) 

1140 IF J=I THEN 1160 

1150 LET XCI,J)=XCI,J)/XCI,I) 

1160 NEXT J 

1170 LET XCI,I)=1 

1180 FOR L=0 TO K 

1190 IF L=I THEN 1270 

1200 LET YCL)=YCL)-XCL,I)*YCI) 

1210 FOR J=0 TO K 

1220 LET ACL,J)=ACL,J)-XCL,I)*ACI,J) 

1230 IF J=I THEN 1250 

1240 LET XCL,J)=XCL>J)-XCL,I)*XCI,J) 

1250 NEXT J 

1260 LET XCL,I)=0 

1270 NEXT L 

1280 NEXT I 

1290 LET S6= CCECM),ECM)) 


17 





1300 FOR 1 = 0 TO K 

1310 LET S6=S6-YCI>*CCECI)*ECM)) 

1320 NEXT I 

1330 LET S7=S6/CN-M> 

1340 LET R2=1-S7/CSCM)*SCM>) 

1350 LET R=S0RCR2) 

1360 LET S8=SQRCS7) 

1370 IF P1=0 THEN 1390 

1380 PRINT “VARIANCE-COVARIANCE MATRIX" 

1390 FOR 1 = 0 TO K 
1400 FOR J = 0 TO K 
1410 LET ACI>J>=AC l> J>*S7 
1420 IF P1=0 THEN 1440 
1430 PRINT ACI*J>* 

1440 NEXT J 

1450 IF P1=0 THEN 1480 

1460 PRINT 

1470 PRINT 

1480 NEXT I 

1490 PRINT 

1500 PRINT “INDEX'S "B'S “STD. ERROR'S "T-RATIO" 

1510 FOR 1 = 0 TO K 

1520 PRINT ECI)> YCI)> SQRCACI,I))> YCI5/SQRCAC I,I>> 
1530 NEXT I 
1540 PRINT 

1550 PRINT “R-SQUARED 3 " R2.» “R=“ R 
1560 PRINT 

1570 PRINT “STAND. ERROR OF EST. =“i S8> “D. F. = "* CN-M) 

1580 PRINT 

1590 FOR 1=1 TO N 

1600 LET Z=DCI>ECM))-YCO> 

1610 FOR J=1 TO K 

1620 LET Z=Z-YCJ)*DCI,ECJ>> 

1630 NEXT J 
1640 LET QCI)=Z 
1650 NEXT I 
1660 LET W=0 
1 670 FOR 1 = 2 TO N 

1680 LET W=W+CQCI)-Q(1-1)>*CQCI)-Q(I-1)) 

1690 NEXT I 
1700 PRINT 

1710 IF P2=0 THEN 1780 

1720 PRINT "ACTUAL'S “PREDICTED'S “RESIDUAL" 

1730 LET 1=0 
1740 LET 1=1+1 

1750 PRINT DC I * ECM))* DC I » ECM))-Q<I )» QCI) 

1760 IF I = N THEN 1780 
1770 GO TO 1740 
1780 PRINT 

1790 PRINT "DURBIN-WATSON STAT.=“ W/S6 

1800 IF H<G THEN 1820 

1810 GO TO 1850 

1820 PRINT 

1830 GO TO 500 

1840 PRINT "CORRELATION MATRIX BECOMING SINGULAR" 
1850 PRINT 

1860 PRINT " *****PR0BLEM COMPLETED*****" 

1870 STOP 
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1940 DATA 9*4*4 

1950 DATA 247 69*34181*3551 6,29566 
1952 DATA 34181*35516*29566*37009 


1954 DATA 35516*29566*37009*34542 
1956 DATA 

295 66* 37009*34542*41485 

37007* 34542* 41485* 49972 
9*4,4 

24769* 34181*35516*29566 
34181* 3551 6* 29566* 37009 
3551 6*29566*37009*34542 
29566,37009* 34542* 41485 
37007* 34542* 41485* 49972 
34542* 41485* 49972* 53944 
41485*49972* 53944* 47062 
49972* 53944* 47062*47284 
53944* 470 62* 47284* 57765 
1 * 1 * 1 * 1 
3*4 

2 * 2 * 1,0 
1*3* 4 
3*3*0*1 
1 *2* 3* 4 
4* 1 * 0* 0 
4*3 


1958 

DATA 

1960 

DATA 

1961 

DATA 

1962 

DATA 

1963 

DATA 

1964 

DATA 

1965 

DATA 

1966 

DATA 

1967 

DATA 

1968 

DATA 

1969 

DATA 

1970 

DATA 

1972 

DATA 

1974 

DATA 

1980 

DATA 

1982 

DATA 

1984 

DATA 

1990 

DATA 

1992 

DATA 

2046 

END 
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1* BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 
2' MODIFIED VERSION OF DARTMOUTH STATISTICS PROGRAM 
3* DARTMOUTH PROGRAM NAME: R-X-C*** 

4* MODIFIED: JULY 2 1970 

5* 

6* B.C.C..LIBRARY PROGRAM NAME: RXC 

7' 

8* DESCRIPTION-- ANALYSIS OF TABLE WITH R ROWS AND C COLUMNS 

9* 

10* SOURCE-- UNKNOWN 

1 1 » 

12* INSTRUCTIONS-- TYPE RUN FOR INSTRUCTIONS. 

13* 

14* 

15* 

1 6* 

70 READ T 
75 DATA 8 
80 DIM 0C6>,HC6) 

85 DIM DC 6, 6),EC36>, F( 6, 6), RC 6)* CC 6), SC 6, 6) 

90 DIM QC 6) , UC 6) 

95 PRINT "HOW MANY ROWS DOES YOUR TABLE HAVE": 

100 INPUT A 

105 IF A> 6 THEN 125 

110 PRINT "COLUMNS": 

115 INPUT B 

120 IF B< = 6 THEN 140 

125 PRINT "THIS PROGRAM WILL NOT ACCOMMODATE A TABLE LARGER THAN" 
130 PRINT "6 X 6. RE-ENTER NUMBER OF ROWS" 

135 GO TO 100 

140 PRINT "PLEASE TYPE IN YOUR DATA." 

144 PRINT "ENTER TABLE BY ROW, ONE ENTRY PER LINE " 

145 FOR I = 1 TO A*B 

146 INPUT EC I) 

147 NEXT I 
150 PRINT 
155 PRINT 
160 PRINT 
165 PRINT 
170 PRINT 

175 LET N = A*B 

130 IF N=CA*B> THEN 200 

185 PRINT "YOU HAVE NOT ENTERED THE DATA CORRECTLY, BUT THAT'S" 
190 PRINT "ALL RIGHT, CHIEF, JUST TRY AGAIN." 

195 GO TO 145 

200 PRINT "YOUR ",*A;"BY "J B5 "TABLE:" 

205 PRINT 

210 FOR I = 1 TO A 
215 FOR J = 1 TO B 

220 LET RCI) = 0 

221 LET CCJ) = 0 
225 NEXT J 

230 NEXT I 
235 LET N = 0 
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240 FOR I = 1 T3 A 
245 F3R J = 1 TO B 
250 LET DCI,J> = ECCI-1>*B + J) 

255 LET R(I) = R(I) +"DCI#J> 

260 LET CCJ) = CCJ) + DC I« J) 

265 LET N = N + DCI..J) 

270 PRINT DCI,J); TABCJ*T>; 

275 NEXT J 

280 PRINT "J "; RCI) 

285 NEXT I 

290 FOR I = 1 TO CCB*T + 13)/2) 

295 PRINT "; 

300 NEXT I 
305 PRINT 

310 FOR J = 1 TO B 

315 PRINT CCJ); TABCJ*T); 

320 NEXT J 
325 PRINT M : M ;N 
330 PRINT 
335 PRINT 

340 PRINT "IS THIS DATA TABLE CORRECT"; 

345 INPUT AS 

350 IF AS = "NO" THEN 370 
355 IF AS = "YES" THEN 380 
360 PRINT "WHAT"; 

365 GO TO 345 

370 PRINT "THEN TYPE THE TABLE IN AGAIN" 

375 GOTO 144 

380 FOR I = 1 TO A 

385 FOR J = 1 TO B 

390 LET FCI,J) = RCI)*CCJ)/N 

395 LET SCI,J) = CDCI,J) - FCI,J))t2/FCI,J) 

400 LET C2 = C2 + SCI>J) 

405 NEXT J 

410 NEXT I 

415 LET G1 = C2 

420 LET M = CA-1)*CB-1) 

425 LET G = Gl/M 

430 LET P = 1 

435 LET N9 = 1000 

440 IF G < 1 THEN 465 

445 LET A9 = M 

450 LET B9 = N9 

455 LET F = G 

460 GO TO 480 

465 LET A9 = N9 

470 LET B9 = M 

475 LET F = 1/G 

480 LET A1 = 2/C9*A9) 

485 LET B1 = 2 /C9*B9) 

490 LET Z = ABSC C1-B1)*F»C.333333)- 1+A1) 

495 LET Z = Z/SQRC Bl*Ft C• 666667)+Al) 

500 IF B < 4 THEN 520 

505 LET P=<1+Z*(.196854+Z*(.115194+Z*C.000344+Z*.019527))))f4 
510 LET P = .5/P 
515 GO TO 530 

520 LET Z = Z*<1 + .08*Zt4/Bt3) 

525 GO TO 505 

530 IF G < 1 THEN 540 

535 GO TO 545 
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540 LET P = 1-P 

545 LET PI = INTC10000%P)/10000 
550 IF PI = 0 THEN 565 
555 PRINT 
560 GO T0 590 

565 PRINT "THE SIGNIFICANCE LEVEL IS LESS THAN .0001. THIS IS " 
570 PRINT "MUCH RARER THAN HEN’S TEETH IN RANDOM NUMBERS. ” 

575 PRINT "IN FACT# YOU MIGHT SAY RARER THAN RUM’S TEETH IN" 

580 PRINT "HANDOM RUMBERS!" 

585 GO TO 715 

590 PRINT "THE SIGNIFICANCE LEVEL IS ”J PI 

595 PRINT 

600 LET P = PI 

605 IF P > .1 THEN 635 

610 IF P > .05 THEN 685 


615 IF P > .01 THEN 695 
620 IF P > .001 THEN 705 

625 PRINT "P < OR = .001. RARER THAN HEN’S TEETH IN RANDOM NUMBERS!" 
630 GO TO 715 

635 PRINT "P > .10. GET MORE CASES OR FORGET IT." 

640 PRINT "D3 YOU WANT MORE OUTPUT <PERCENTAGE TABLES# GAMMA# ETC. >"# 
645 INPUT AS 

650 IF AS = "NO" THEN 1640 
655 IF AS = "YES" THEN 670 
660 PRINT "WHAT"# 

665 GO TO 645 


670 PRINT "O.K. HERE GOES# 
675 GO TO 715 
680 GO TO 715 


685 PRINT 
690 GO TO 
695 PRINT 
700 GO TO 
705 PRINT 
710 PRINT 
715 PRINT 


".10 >= P > .05 
71 5 

".05 >= P > .01. 
715 

".01 .= P > .001 

"RANDOM NUMBERS. 


it 


BUT I WOULDN’T GIVE YOU TWO CENTS FOR IT.". 

BORDERLINE. ONE IN 11 ISN'T ALL THAT RARE 
MOST PEOPLE WOULD BUY 20 TO 1 ODDS." 

YOU DON *’T SEE MANY LIKE THAT WITH " 


720 PRINT 

725 PRINT "WOULD YOU LIKE TO KNOW THE ABSOLUTE VALUE OF CHI-SQUARE"; 
730 INPUT AS 

735 IF AS = "YES" THEN 755 

740 IF AS = "NO" THEN 775 

745 PRINT "COME ON# YES OR NO"; 

750 GO TO 730 


755 PRINT 
760 GOSUB 
7 65 PRINT 
770 PRINT 
775 PRINT 
780 PRINT 
785 PRINT 
790 INPUT 
795 IF AS 
800 PRINT 
805 PRINT 
810 INPUT 
815 IF AS 
820 IF AS 
825 IF AS 
835 PRINT 
840 GO TO 


1 665 

TABC 9 ) # 

"CHI-SQUARE IS "#ABSC C2) ;"WITH (A-1)*(B-1);"D.F." 


"DO YOU WANT THE PERCENTAGE TABLE"; 

AS 

= "NO" THEN 1065 

"TYPE ’H’ IF YOU WANT IT DONE HORIZONTALLY# ’V’ IF" 
"VERTICALLY# AND ’B’ IF BOTH." 

AS 

= "H" THEN 940 
= "V" THEN 845 
= "B" THEN 8 45 
"WHAT"; 

810 
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845 PRINT " VERTICAL PERCENTAGE TABLE" 

850 PRINT 

855 F3R I = 1 T0 A 
8 60 FOR J = 1 TO B 

865 LET U1 = INTC1000*CDCI,J)/CCJ>>+.5)/10 
870 LET UCJ) = UCJ) + U1 
875 PRINT Ul;TABCJ*T>j 
880 NEXT J 

885 PRINT "J INT(1000*RC I ) /N + • 5 )/1 0 

890 LET V = V + INTC1000*RCI)/N + .5)/10 
895 NEXT I 

900 FOR I = 1 TO CCB*T + 13)/2) 

905 PRINT 
910 NEXT I 
915 PRINT 

920 FOR J = 1 TO B 

925 PRINT UCJ>* "%'S TABCJ*T)i 

930 NEXT J 

935 GO TO 945 

940 GO TO 9 60 

945 PRINT "5V5"%" 

950 PRINT 
955 I 

F A$ = "V" THEN 1065 

960 PRINT " HORIZONTAL PERCENTAGE TABLE'S* 

965 PRINT 

970 FOR I = 1 TO A 
975 FOR J = 1 TO B 

930 LET T1 = INT(1000*CDCI»J)/RCI>)+.5>/10 
985 LET QCI)=Q(I) + T1 
990 PRINT Tl;TABCJ*T); 

995 NEXT J 
1000 PRINT 
1005 NEXT I 

1010 FOR I = 1 TO CCB*T + 13>/2> 

1015 PRINT 

1020 NEXT I 

1025 PRINT 

1030 FOR J = 1 TO B 

1035 LET H5=INTC1000+(CCJ)/N)+.5)/10 
1040 LET 176 = W6 + H5 
1045 PRINT H5; TABCJ*T); 

1050 NEXT J 

1055 PRINT "J "} W65 

1060 PRINT 

1065 PRINT 

1070 PRINT 

1075 PRINT 

1080 PRINT "DO YOU WANT THE EXPECTED CELL VALUES' 
1035 INPUT AS 

1090 IF AS = "NO" THEN 1220 
1095 IF AS = "YES" THEN 1110 
1100 PRINT"WHAT"i 
1105 GO TO 1035 
1110 PRINT 

1115 PRINT TABC5) i "EXPECTED CELL VALUES" 

1120 PRINT 
1125 PRINT 
1130 FOR I = 1 TO A 
1135 FOR J = 1 TO B 

1140 PRINT INTC100*FCI*J)+.5)/100JTABCJ+T)5 
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1145 LET HCI) = HCI) + INK 1 00* FC I > J ) +.5)/100 
1150 LET 0 C J ) = 0CJ) + INK100*FC l, J) +.5)/100 
1155 NEXT J 

1160 PRINT ”5 HCI) 

1165 NEXT I 

1170 F0R I = 1 T3 CCB*T + 13>/2> 

1175 PRINT ". "; 


1180 NEXT I 

1185 PRINT 

1190 FOR J = 1 TO 3 

1195 PRINT OCJ); TABCJ*T); 

1200 NEXT J 

1205 PRINT "J "; N 

1210 PRINT 

1215 PRINT 

1220 PRINT 

1225 PRINT 

1230 PRINT M D0 YOU WANT THE DIFFERENCE BETWEEN THE OBSERVED VALUE" 
1235 PRINT "AND THE EXPECTED VALUE FOR EACH CELL"; 

1240 INPUT AS 

1245 IF AS = "YES" THEN 1265 
1250 IF AS = "NO" THEN 1355 
1255 PRINT "WHAT"; 

1260 GO TO 1240 
1265 PRINT 
1270 PRINT 


1275 
1280 
1285 
1290 
1295 
1300 
1305 
1310 
1315 
1320 
1325 
1330 
1335 
1340 
1345 
1350 
1355 
1360 
1365 
1370 
1375 
1380 
1385 
1390 
1395 
1400 
1405 
1410 
141 5 
1420 
1425 
1430 
1435 
1440 


PRINT 

PRINT 

PRINT " OBSERVED VALUE MINUS EXPECTED VALUE" 

PRINT 

FOR I = 1 TO A 
FOR J = 1 TO B 

PRINT INK 100*CDCI>J)-FCI>J)))/100; TABCJ*T); 

NEXT J 
PRINT ": " 

NEXT I 

FOR I = 1 TO CCB*T + 7)/2) 

PRINT". 

NEXT I 
PRINT 

PRINT TABC B*T);":" 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT "ARE BOTH THE ROW AND COLUMN SCALES ORDINAL"; 
INPUT AS 

IF AS = "NO" THEN 1650 
IF AS = "YES" THEN 1405 
PRINT "WHAT"; 

GO TO 1380 

PRINT "DO YOU WANT TO KNOW GAMMA"; 

INPUT AS 

IF AS = "NO" THEN 1595 
IF AS = "YES" THEN 1440 
PRINT "WHAT DID YOU SAY" 

PRINT "PLEASE TYPE 'YES’ OR ’NO’" 

GO TO 1410 
LET G1 = 0 
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1445 FOR I = 1 TO CA - 1) 

1450 FOR J = 1 TO (B - 1) 

1455 LET X = 0 
1460 FOR K = (1 + 1 > TO A 
1465 FOR L = (J+l > TO B 
1470 LET X = X + DCi<*L) 

1475 NEXT L 
1480 NEXT K 

1485 LET G1 = G1 + <D<I«J)*X> 

1490 NEXT J 

1495 NEXT I 

1500 LET G2 = 0 

1505 FOR I = 1 TO <A-1 ) 

1510 FOR J = B TO 2 STEP -1 

1515 LET X = 0 

1520 FOR K = <1 + 1 ) TO A 

1525 FOR L = <J-1> TO 1 STEP -1 

1530 LET X = X + D(K*L> 

1535 NEXT L 
1540 NEXT K 

1545 LET G2 = G2 + CD<I*J)*X> 

1550 NEXT J 
1555 NEXT I 

1560 LET G9 = CG1-G2) / (G1+G2) 

1565 PRINT 
1570 GOSUB 1665 

1575 PRINT TABC20);"GAMMA IS "J -INTC10000*G9)/10000 
1530 PRINT 
1585 PRINT 

1595 PRINT ^TsS-TSK. NO INTELLECTUAL CURIOSITY." 

1600 PRINT "ARE BOTH THE ROWS AND THE COLUMNS INTERVAL SCALES"; 

1605 INPUT AS 

1610 IF AS = "NO" THEN 16 60 
1615 IF AS = "YES" THEN 1630 

1620 PRINT "HOW’S THAT AGAIN? A SIMPLE ’YES’ OR ’NO’ WILL SUFFICE. 
1625 GO TO 1605 

1630 PRINT "HAVE YOU CONSIDERED PRODUCT MOMENT CORRELATIONS? I" 
1635 PRINT "CAN'T DO THEM* BUT THERE ARE OTHER PROGRAMS WHICH CAN." 
1640 PRINT "SO LONG. 

1645 STOP 

1650 PRINT "THAT’S ALL I CAN DO FOR YOU. BYE-BYE." 

1655 STOP 

1660 PRINT "THEY AREN’T? THEN I'LL BE SHOVING ALONG. GOOD-BY." 

1665 FOR I = 0 TO 9 

1670 PRINT " "5 

1675 NEXT I 

1635 REM 

1690 REM 

1695 RESTORE 

1700 REM 

1705 END 

READY 
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10 ’ BERKSHIRE COMMUNITY COLLEGE USER PROGRAM LIBRARY 

20 * MODIFIED VERSION OF DARTMOUTH COLLEGE STATISTICS PROGRAM 

30 * DARTMOUTH PROGRAM NAMES CORREL*** 

40 * MODIFIED: JULY 2 1970 
90 ’ 

100 ' B.C.C. LIBRARY PROGRAM NAME: C0RREL 
110* 

120' DESCRIPTION--COMPUTES A CORRELATION COEFFICIENT 
130’ 

140* SOURCE—UNKNOWN 

150' 

1 60 * INSTRUCTIONS--THE DATA CONSISTS OF N SETS OF NUMBERS. 

170* EACH SET CONSISTS OF TWO ELEMENTS* X AND Y. N IS 
180* ENTERED AS DATA IN LINE 900. IN LINE 910 AND FOLLOWING 
190* ENTER DATA ALTERNATING X AND Y VALUES. SAMPLE DATA 
200* ARE IN LINES 900-910. 

210' 

220' ****** MAIN PROGRAM * * ******* 

230’ 

240* 

250 READ N 
260 FOR I = 1 TO N 
270 READ X* Y 
280 LET SI = SI + X 

290 LET S2 = S2 + Y 

300 LET S3 = S3 + X*Y 

310 LET S4 = S4 + X*X 

32u LET S5 = S5 + Y*Y 

330 NEXT I 

340 PRINT "THE CORRELATION COEFFICIENT = "i 
350 LET A = N*S3 - S1*S2 

360 LET B = SQRC(N*S4 - S1»2)*CN*S5 - S2*2)) 

370 LET A = INTCA/B*1000+.5)/1000 
380 PRINT A 
900 DATA 5 

910 DATA 1 * 5* 2* 3* 3* 0* 4* - 5* 5* - 1 1 
2046 END 

READY 
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